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Abstract 

Involvements of major financial institutions in the recent financial crisis have generated renewed 
interests in fragilities of global financial networks among economists and regulatory authorities. In 
particular, one potential vulnerability of the financial networks is the "financial contagion" process in 
which insolvencies of individual entities propagate through the "web of dependencies" to affect the entire 
system. In this paper, we formalize a banking network model originally proposed by researchers from 
Bank of England and elsewhere that may be applicable to scenarios such as the OTC derivatives market, 
define a global stability measure for this model, and comprehensively evaluate the stability measure 
over more than 700,000 combinations of networks types and parameter combinations. Based on such 
comprehensive evaluations, we discuss some interesting implications of our evaluations of this stability 
measure, and derive topological properties and parameters combinations that may be used to flag the 
network as a possible fragile network. 
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1 Introduction 



The recent global financial crisis at an unprecedented level have exposed potential weaknesses of the global 
economic system, renewing interests in the determination of fragilities of various segments of the global econ- 
omy. Since financial institutions played a crucial role in this crisis, financial system governed by borrowing, 
lending and participation in risky investments have attracted a major part of the attention of economists; 
see [1, 2] for a survey. The issue of instability of free market based financial systems is not new and has been 
under discussion among the economists starting with the early works of Fisher [3] and Keynes [4] during 
the 1930's great depression era. However, the exact causes of such instabilities have not been unanimously 
agreed upon yet. For example, economists Ekelund and Thornton contend that a major reason for the recent 
financial crisis is the enactment of the Gramm-Leach-Bliley Act in 1999 that removed several restrictions on 
mixing investment and consumer banking [5], whereas many other economists such as Calabria disagree with 
such an assertion [6]. Some economists such as Minsky in fact have argued that such instabilities are are 
systemic (i. e., inherent) for many modern market-based economic systems [7]. 

Our motivation in this paper to investigate the issue of global stabilities of financial networks comes from 
the point of view of a regulatory agency (as was also done, for example, in [1]). A regulatory agency with 
sufficient knowledge about a significant part of a global financial network is expected to periodically evaluate 
the stability of the network, and flag the network ex ante for further analysis if it fails some preliminary test 
or exceeds some minimum threshold of vulnerability. In this motivation, flagging a network as vulnerable 
does not necessarily imply that such is the case, but that such a network requires further analysis based on 
other aspects of free market economics that are not or simply cannot be modeled 1 . While too many false 
positives may drain the finite resources of a regulatory agency for further analysis and investigation, this 
motivation assumes that vulnerability is too important an issue to be left for an ex post analysis. 

2 Informal overview of our contribution 

Although simpler topological properties such as clustering coefficients have been used by authors to study 
properties of financial network [8,9], they are too simplistic to be useful in the type of stability analysis done 
in this paper. Here we formalize a simplified ex ante financial network contagion model similar to what has 
been used by prior researchers from Bank of England and elsewhere [1, 10-12], and define a global stability 
measure for this model. We perform a thorough evaluation and analysis of this stability measure over more 
than 700,000 combinations of networks types and parameters. Such a comprehensive analysis allows us to 
uncover many interesting insights into the relationships of the stability with other relevant parameters of 
the network, such as: 

Effect of unequal asset distribution: More uneven distribution of assets among banks in a banking 
network tend to make the network more vulnerable. Moreover, a banking network with uneven 
distribution of assets among banks tends to exhibit the following properties: 

• Failures of those banks with disproportionately higher assets contribute more damage to the 
long-term stability of the network. 

1 For example, some such factors are the rumors and panics caused by the insolvency of a large bank and a possible subsequent 
credit freeze. While fears, panics and rumors are all real aspects in networked economics, there are hardly any universally agreed 
upon good way of modelling them. 
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• These networks have a minimal instability even if their equity to asset ratio is large and 
comparable to loss of external assets. This is not the case for networks where the distribution 
of assets among banks is uniform. 

Thus, in summary, we conclude that banks with disproportionately large external assets affect 
the stability of the entire network in an adverse manner. 

Effect of total amount of assets: Stabilities of networks with highly uneven distribution of assets 
among banks are affected very little by the amount of total external assets. 

Effect of connectivity: For banking networks where all banks have the same amount of external 
assets, higher connectivity leads to lower stability. In contrast, for banking networks in which few 
banks have disproportionately higher external assets compared to the remaining banks, higher 
connectivity leads to higher stability. 

Correlated versus idiosyncratic insolvencies: Correlated initial insolvencies of banks cause sig- 
nificantly more damage to the stability of a banking network as opposed to idiosyncratic initial 
insolvencies. 

Phase transition properties of stability: The stability measure exhibits several sharp phase tran- 
sitions for various banking networks and combinations of parameters. We have also provided 
intuitive theoretical explanation for one such phase transition. 

Although the issue of stability of financial systems has been discussed by prior researchers [1, 10-17, 19, 20, 23], 
to our knowledge no prior paper has mathematically defined a global stability measure and performed a 
comprehensive evaluation of such a measure as is done in this paper. 

2.1 Policy implications 

Returning to our original motivation of flagging financial networks for potential vulnerabilities, our empirical 
results suggest that a network model similar to that used in the paper may be flagged for the following cases: 

• the equity to asset ratio of banks in the network is extremely low, or 

• the network has a highly skewed distribution of external assets and inter-bank exposures among its 
banks and the network is sufficiently sparse, or 

• the network does not have either a highly skewed distribution of external assets or a highly skewed 
distribution of inter-bank exposures among its banks, but the network is sufficiently dense. 

The reader should note that the above policy implications are valid only in the context of the specific network 
model and parameter ranges considered in the paper. 

3 Network model and global stability measure 

An investigation of the stability of a financial system may begin by addressing the following items: 

network model: How do we precisely model the individual components and their mutual interactions 
in the financial system under consideration ? 

insolvency propagation: How does insolvencies of individual components of the financial system 
affect the overall network over time ? 
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global stability evaluation: How do we evaluate the global stability of the financial system ? 
We describe these items in the sequel. 

3.1 The network model 

As was done by prior researchers, for simplicity we formulate our results in terms of balance-sheet insolvency 
cascades in a network of financial institutions (hereafter called "banks" and "banking networks") with 
interlinked balance sheets, where losses flow into the asset side of the balance sheets. The same formulation 
can be reused with appropriate modifications for analyzing cascades of cash-flow insolvency in OTC markets. 
From now on, we refer to balance-sheet insolvency simply as insolvency. 

A banking network is represented by a node- weighted and edge- weighted directed graph G = (V, F) in 
which the node set V represents the banks and the edge set F represents direct inter-bank exposures. Let 
the notation in-degree(w) denote the in-degree of node v £ V, and the notation w(e) = w(u, v) > denote 
the weight of an edge e = (u, v) G F. Throughout, we use E and I to denote the total external asset and the 
total inter-bank exposure of the banking network, respectively. Two types of banking networks are possible: 

homogeneous network model: E and I are equally distributed between the nodes and the links, 
respectively, and 

heterogeneous network model: E and / may not be equally distributed between the nodes and 
the links, respectively. 

For both types of networks, we define the following quantities: 

• a v denotes the share of the total external asset E for bank v eV. 

• 7 G (0, 1) denotes the ratio of equity to asset. 

• t'v = J2( v u)eF w ( v > u ) denotes the total interbank asset of bank v G V. 

• K — J2( u v)ef w ( u , v ) denotes the total interbank borrowing of bank ueV. 

• e v = (b v — l v ) + cr v (E — J2vev — Lv )) denotes the effective share of the total external asset E for 
bank v G V. 

• a v = e v + l v denotes the total asset of bank v G V. 

• c v — 7 a v denotes the net worth (equity) of bank v G V. 

Both homogeneous and heterogeneous network models are relevant in practice, and have been investigated 
by prior researchers. Calculations of these parameters for a simple 5-node banking network are shown in 
Section A in the appendix. 

3.2 Initial insolvency via shocks 

The initial insolvencies of a banking network at time t — are caused by "shocks" received by a subset of 
banks. Such shocks can occur, for example, due to operational risks (e. g., frauds 2 ) or credit risks, and has 
the effect of reducing the external assets of a selected subset of banks. Mathematically, we are given the 

2 Studies such as by Iyer and Peydro [34] show that fraud is an important cause of bank insolvency. 
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real number $ e (0, 1] > 7 denoting the severity of the initial shock, and the effect of the initial shock is to 
simultaneously decrease the external assets of each shocked bank v from e v by s v = $ e v , thereby reducing 
the net worth of v from its original value c v to its new value c v — s v . By the phrase "shocking mechanism", 
we refer to the rule that is used to select the initial subset of nodes to be shocked. 

3.3 Insolvency propagation equation 

The insolvencies are assumed to propagate in discrete time units, say t = 0, 1,2, We add "(i)" to all 

relevant variables to indicate their dependences on time. A bank becomes insolvent if its modified net worth 
becomes negative, and such a bank is removed from the network in the next time step. Let insolvent (t) C V 
be the set of banks that have become insolvent before time t. The insolvency of a bank at time t affects the 
equity of other banks in the network by the following discrete-time "insolvency propagation equation" : 



mini I c v (t) I , b v \ 

VueV\V insolvcnt (t): c u (t + 1) = c u (t) - V — -. - } 

^— ' m-degreefw, t) 

(c„(t)<0) 

V. A («€V\Vinsolvent(t)) 



(1) 



In the above equation, the term : j Ct, ^j allows the loss of equity of an insolvent bank to be distributed 

in — degree(v, t) 

b v 

equitably among its creditors that have not become insolvent yet, whereas the term ; r- ensures 

in — degree(v, t) 

that the total loss propagated cannot be more than the total interbank exposure of the insolvent bank. The 
Insolvency propagation continues until no new bank becomes insolvent. A pseudo-code of the insolvency 
propagation method starting with the initial shock is shown in Fig. 1. 

3.4 Rationale for the network and insolvency propagation model 

As several prior researchers such as [10,23,24] have commented: 

"conceptual frameworks from the theory of weighted graphs with additional parameters may provide 
a powerful tool for analysis of banking network models". 

Several parametric graph-theoretic models have been used by prior researchers in finance, economics and 
banking industry to study various properties and research questions involving banking systems [1, 10, 19, 20, 
25-30], differing in the way edges are interpreted and additional parameters are used to characterize the 
contagion. As noted by researchers such as [10,29]: 

"the modelling challenge in studying banking networks lies not so much in analyzing a model that is 
flexible enough to represent all types of insolvency cascades, but in studying a model that can mimic 
the empirical properties of these different types of networks". 

The insolvency propagation model formalized and evaluated in this paper using a mathematically precise 
abstraction is similar in spirit to the models used in [1, 10-12, 25], and is represented by cascades of cash-flow 
insolvencies. As observed in [29], OTC derivatives markets and similar other markets are prone to such type 
of cascades. In such markets parties deal directly with one another rather than passing through an exchange, 
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and thus each party is subject to the risk that the other party does not fulfill its payment obligations. The 
following example from [29] illustrates chains of such interactions: 

[29] "Consider two parties A and B, such that A has a receivable from party B upon the realization 
of some event. If B does not dispose of enough liquid reserves, it will default on the payment. Now 
consider that B has entered an off-setting contract with another party C, hedging its exposure to the 
random event. If C is cash-flow solvent, then the payment will flow through the intermediary B and 
reach A. However, if C is cash-flow insolvent and defaults, then the intermediary B might become 
cash-flow insolvent if it depends on receivables from C to meet its payment obligations to A" . 

The length of such chains of interactions in some OTC markets, like the credit default swap market, is 
significant [21, 22], thereby increasing the probability of a cascade of cash- flow insolvencies [23]. As observed 
in [10] and elsewhere, a insolvency propagation model such as studied in this paper 

"conceptualises the main characteristics of a financial system using network theory by relating the 
cascading behavior of financial networks both to the local properties of the nodes and to the underlying 
topology of the network, allowing us to vary continuously the key parameters of the network" . 

Although the cascading effect studied is a simplified one, as noted by Haldane and May from Bank of England 
in their paper [1]: 

"This is a deliberate oversimplification, aimed at a clearer understanding of how an initial failure can 
propagate shocks throughout the system" . 

3.5 The global stability measure 

Consider a banking network with all parameters as described before, and let < K, < 1 be a real number 
denoting the fraction of nodes in V that received the initial shock under a shocking mechanism T. The 
vulnerability index of the network is then defined as 

£ (£> G, 7, <&, T) = -j— x max ( lim I insolvent (i) I } 

|V| V'CV: \V'\=K L t^co 1 1 J 

For example, £ (0.1, G, 0.3, 0.5, T) = 0.9 means it is possible to make 90% nodes of the network G become 
insolvent with 7 = 0.3 and $ = 0.5 if we provide an initial shock to a suitably selected subset of 10% of 
nodes of G under the shocking mechanism T. Note that lower values of £ imply higher global stability of a 
network. For simplicity, we will omit the arguments of £ when they are clear from the context. 

3.6 Rationale for the global stability measure 

It is possible to think of other alternate measures of global stability than the one quantified above. However, 
the measure introduced above is in tune with the ideas that references [1, 10,25] directly (and, some other 
references such as [12, 26, 27] implicitly) used to empirically study their networks. Thus, we decided to follow 
the cue provided by other researchers in the banking industry in studying various insolvency propagation 
models in defining our stability measure. Similar kinds of measures have also been used by prior researchers 
in social networks in the context of influence maximization issues [31, 32]. 
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3.7 Network topology generation 

We consider two topology models previously used by economists to generate random banking networks: 



• the directed scale- free (SF) network [33] that has been used by prior banking network researchers such 
as [29,30,35,36], and 

• the directed Erdos-Renyi (ER) network [41] that has been used by prior banking network researchers 
such as [11,37-40]. 

The directed scale-free networks in this paper are generated using the algorithm outlined by Bollobas et 
al. [41]. The graph grows by adding single edges at discrete time steps. Let o, /3, rj 7 Si n (in-degree) and S ou t 
(out-degree) be non- negative real numbers with a + (3 + r\ = 1. The initial graph G(0) at time t = has 
just one node with no edges. At any time t the graph G(t) has exactly t edges, and a random number n(t) 
of nodes. For t > 0, the graph G(t + 1) is obtained from the graph G(t) by using the following rules: 

• With probability a, add a new node v together with an edge from v to an existing node w, where w 
is chosen randomly such that Pr{w = u] = ^t+i^wft)")^ ^ or everv existing node u, where di n (u) is the 
in-degree of node u in G(t). 

• With probability (3, add an edge from an existing node v to an existing node w, where v and w are 
chosen independently, such that 



with d out (u) being the out-degree of node u in G(t). 
• With probability 77, add a new node w and an edge from an existing node v to w, where v is chosen 



To study the affect of connectivity on network stability, we generated random SF and ER networks with 
average degrees of 3 and 6. In addition, to study the effect of sparse hierarchical topology of networks 
on its stability, we used the Barabasi-Albert preferential- attachment SF model [33] to generate random in- 
arborescence networks. In-arborescences are directed rooted trees with all edges oriented towards the root, 
and have the following well-known topological properties: 

• They belong to the class of sparsest connected directed acyclic graphs. 

• They are hierarchical networks, i. e., the nodes can be partitioned into levels L\, L2, . ■ . , L p such that 
L\ has exactly one node (the "root") and nodes in any level Li have directed edges only nodes in level 

(see Fig. 2). In particular, the root models a "central bank" that only lends money to other banks 
but does not borrow money from any bank. 

A random in-arborescence network of n nodes is generated using Barabasi-Albert preferential-attachment 
model [33] in the following manner: 

• Start with one node as the root of the tree. 



Pr[v = u] 
Pr[w = u] 




- for every existing node u 
for every existing node u 
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• At every successive time step, add a new node to the tree, and add an undirected edge from the 
new node to a randomly selected node u that is already in the tree. Node u is selected such that 
Pr[u — v] — ^-■ dC jcgrce( M ,) i where degree(y) denotes the degree of node y. Increment the degree of node 
u and the new node. Repeat this step till the tree has n nodes. 

• Orient all the edges towards the root. 

The minimum difference between two non-identical values of the average vulnerability index over 10 runs for 
two n-node networks is 1 /(Wn). Thus, to allow for minor statistical biases introduced by any random graph 
generation method, we consider two vulnerability indices to be same (within the margin of statistical error) 
if their absolute difference is no more than 1 /(3n), which is above 1 /(Wn) but no more than 0.7% of the total 
range of the vulnerability indices. For example, if £i and £2 are the average vulnerability indices (over 10 
runs) of two 50-node networks then we consider £1 to be at least £2 if £1 > £2 — 0.0066, and consider £1 to 
be the same as £2 if | £1 — £2) < 0.0066. 

3.8 Shocking mechanisms T 

We used the following two mechanisms to select the nodes to receive the initial shock: 

Idiosyncratic (random) shocks: We select a subset of nodes uniformly at random. This corre- 
sponds to the case of random idiosyncratic initial insolvencies of banks, and is a choice that has 
been used by prior researchers such as [10, 25]. 

Coordinated shocks: In this type of non-idiosyncratic correlated shocking mechanism, we seek to 
play an adversarial role in selecting a subset of nodes for the initial shock that may contribute 
more damage to the stability of the network 3 . 

For homogeneous networks, all banks have the same share of the total external asset E. However, 
note the total interbank exposure b v of a bank v is directly proportional to the in-degree of the 
corresponding node v, and, as Equation (1) also suggests, nodes with higher inter-bank exposures 
are the banks that are more likely to transmit the shock to a larger number of other banks. 
Thus, we play an adversarial role by selecting a set of k\V\ nodes in non-increasing order of their 
in-degrees starting from a node with the highest in-degree. 

For heterogeneous networks, nodes with higher "weighted" in-degrees (i. e., higher values of the 
sum of the weights of incoming edges) represent banks that are larger than other banks in terms of 
their external assets, and have more inter-bank exposures. Thus, for heterogeneous networks, we 
play an adversarial role by selecting k\V\ nodes in non-increasing order of their weighted in-degrees 
starting from a node with the highest weighted in-degree. 

The coordinated shocking mechanism falls under the general category of non-idiosyncratic correlated shocks 
where the nodes with high in-degrees or high weighted in-degrees are correlated. 

3.9 Parameter spaces 
3.9.1 Homogeneous networks 

For homogeneous networks, all combinations of the following range of parameters were explored exhaustively: 

3 See, for example [42], about the role of adversarial strategies in measuring the worst-case bounds for network properties. 
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• \V\ = 50,100,300 ; 

• £// = 0.25, 0.5, 0.75, 1, 1.25, 1.5, 1.75, 2, 2.25, 2.5, 2, 75, 3, 3.25, 3.5 ; 

• $ = 0.5,0.6,0.7,0.8,0.9 ; 

• 7 = 0.05, 0.1, 0.15, 0.05 ; and 

• JC = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9. 

Thus, we explored a total of 24,570 combinations of parameters for each of the 2 types of shocks (idiosyncratic 
or coordinated) for each of the 5 network models (in-arborescence, ER of average degree 3, ER of average 
degree 6, SF of average degree 3 or SF of average degree 6), resulting in 245,700 different homogeneous 
banking networks. To correct statistical biases, for each such combinations of parameters, shock types and 
network types, we generated 10 corresponding random networks and computed the average value of the 
vulnerability index over these 10 runs. For generation of ER and SF random networks, we selected the values 
of network generation parameters such that the expected number of edges of the network is 3\V\ or 6|V| 
based on whether we require the average degree of the network to be 3 or 6, respectively. 

3.9.2 Heterogeneous networks 

Recall that in a heterogeneous network it is possible to have a few banks whose external assets or interbank 
exposures are significantly larger than the rest of the banks in the networks, i. e., it is possible to have a 
few banks that are "too big" , and thus heterogeneous networks permit investigation of the effect of such big 
banks on the global stability of the entire network. 

To this end, we define a (a, (3) -heterogeneous network G — {V,F) as one in which we select a random 
subset V' of a\V\ nodes and distribute (3E part of the total external asset E equally among these nodes in 
V'. Let F 1 be the edges involving these a\V\ nodes (i. e., edges in F' have at least one end-point from V'). 
Then we distribute (31 part of the total interbank exposure / equally among a random subset of a\F'\ of 
edges from the edges in F' . The remaining parts of E and /, namely (1 — (i)E and (1 — (3)1 respectively, are 
distributed equally among the remaining (1 — a)\V\ nodes and remaining \F\ — (3\F'\ edges, respectively. 

We performed our simulations for (a, /3)-heterogeneous networks for two pairs of a and (3, namely 
(a, (3) = (0.1,0.95) and (a, (3) = (0.2,0.6), thus we obtained a total of 491,400 different heterogeneous 
banking networks. The combination (a, (3) = (0.1,0.95) corresponds to the extreme situation in which 95% 
of the assets and exposures involve 10% of banks, thus creating a minority of banks that are significantly 
larger than the remaining banks. The other combination (a, (3) — (0.2,0.6) corresponds to a less extreme 
situation in which there are a larger number of moderately large banks. 

Relevant data of all networks used in this paper will be made publicly available in the final version of 
the paper. 

4 Results and discussion 

4.1 Effect of unequal distribution of assets 
4.1.1 On global stability 

We comprehensively compared the values of the vulnerability index £ for homogeneous and heterogeneous 
networks of same types and same values of common global parameters. For this purpose, for the same value 
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of the common parameters |V|, E, n, $ and 7, for the same for network type (ER, SF or in-arborescence) of 
same average degree (6, 3 or 1) and for the same shocking mechanism T (coordinated or idiosyncratic), we 
compared the value of £ for the homogeneous network with the corresponding values of £ for both (0.1, 0.95)- 
heterogeneous and (0.2, 0.6)-heterogeneous networks. The comparison results are tabulated in Table 1, where 
most of the entries are significantly higher than 50% (and in fact, are at least 90%). Thus, we conclude: 

(A) networks where all banks have roughly the same external assets display higher global 
stability over similar networks in which fewer banks have a disproportionately higher 
external assets compared to the remaining banks. 

4.1.2 On residual instability 

For homogeneous networks, if the equity to asset ratio 7 is close enough to the severity of the shock $ 
then the network almost always tends to be perfectly stable, as one would intuitively expect. However, 
the above property is not true in general for highly heterogeneous networks in the sense that, even when 
7 is close to <fr, these networks (irrespective of their topologies and densities) have a minimum amount of 
global instability 4 . In Tables 4- 13 we tabulate residual instabilities for different types of homogeneous and 
heterogeneous networks under coordinated and idiosyncratic shocks. The numbers in these tables show, for 
each combination of network types, \V |, shocking mechanism and values of $ and 7 such that |3> — 7I = 0.05, 
the percentage of networks with this combination for which the vulnerability index £ was less than 0.05, 0.1 
or 0.2. As the reader will observe, all the numbers for heterogeneous networks are significantly lower than 
their homogeneous counter-parts. Thus, we conclude: 

(B) a heterogeneous network of any type and any density, in contrast to its correspond- 
ing homogeneous version, has a non-trivial minimum global instability even if its equity 
to asset ratio is very large and close to the severity of the shock. 

4.2 Effect of total external assets 

The Glass-Steagall Act of 1933 was enacted in USA to control speculation by banks by separating consumer 
banking and investment banking. Later, the Gramm-Leach-Bliley Act of 1999 effectively removed this sep- 
aration between investment and consumer banking, and subsequently in 2007 a collapse of major financial 
institutions was on the horizon. A controversial belief correlates these events by asserting that a repeal 
of Glass-Steagall act allowed a ripple effect of insolvencies of individual banks to many other banks in the 
network. In our setting, the quantity E /i controls the total (normalized) amount of external investments of 
all banks in the network. Thus, varying the ratio E /i allows us to investigate the role of the magnitude of 
total external investments on the stability of our banking network 5 . Our results are tabulated in Table 2. 
All the entries in Table 2 are vanishingly close to 0, showing that heterogeneous networks exhibited very 
small average changes in the vulnerability index £ when E /i was varied. Thus, we conclude: 

4 For graphical illustrations to this phenomena, the reader is referred to supplemental Fig. SI and Fig. S2— S3. For example, 
in supplemental Fig. SI, when 7 is 45% and $ is only 5% more than 7, the vulnerability index £ is approximately for all the 9 
combinations of parameters, but in supplemental Fig. S2— S3 all the 18 networks have a value of £ > 0.1 even when the equity 
to asset ratio 7 is 45% and the severity of the shock is only 5% more than 7. 

5 Visual illustrations of Conclusion (C) are shown in supplemental Fig. S4 for homogeneous networks and in supplemental 
Fig. S5 for heterogeneous networks. 
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(C) for all heterogeneous banking networks under consideration, global stabilities are 
affected very little by the amount of total external assets in the system 6 . 

4.3 Effect of network connectivity 

Although it is clear that connectivity and other similar topological properties of a banking network should 
have a crucial effect on its stability, prior researchers have drawn mixed conclusions on this. For example, 
Allen and Gale [23] concluded that networks with less connectivity are more prone to contagion than networks 
with higher connectivity due to improved resilience of banking network topologies with higher connectivity 
via transfer of proportion of the losses in one bank's portfolio to more banks through interbank agreements. 
On the other hand, Babus [17] observed that, when the network connectivity is higher, liquidity can be 
redistributed in the system to make the risk of contagion lower, and Gai and Kapadia [11] observed that 
higher connectivity among banks leads to more contagion effect during a crisis. The mixed conclusions arise 
because links between banks have conceptually two conflicting effects on contagion, namely, 

• more interbank links may increase the opportunity for spreading insolvencies to other banks, but 

• more interbank links may also provide banks with a type of co-insurance against fluctuating liquidity 
flows. 

As our findings show below, these two conflicting effects have different strengths in homogeneous versus 
highly heterogeneous networks, thus justifying the mixed conclusions of past researchers. 

Recall that in a homogeneous network all banks have the same external asset. Table 14 shows sparser 
homogeneous networks with lower average degrees to be more stable for same values of other parameters. 
Thus, we conclude: 

(D) for homogeneous networks, higher connectivity leads to lower stability. 

In a heterogeneous network, there are banks that are "too big" in the sense that these banks have much 
larger external assets and inter-bank exposures compared to the remaining banks. Table 3 shows that 
for heterogeneous network models denser networks with higher average degree are more stable than the 
corresponding sparser networks of lower average degree for same values of other parameters, specially when 
the heterogeneity of the network is larger (i. e., when a = 0.1 and (3 = 0.95). Thus, we conclude: 

(E) for highly heterogeneous networks, higher connectivity leads to higher stability. 

4.4 Correlated versus idiosyncratic shocks 

For most parameter combinations, our results showed that coordinated shocks, which is a type of corre- 
lated shock, resulted in a insolvencies of higher number of nodes as opposed to idiosyncratic shocks for the 
same network with the same parameters, often by a factor of two or more. Statistics of relevant compar- 
isons for various network type are shown in Table 15. For example, Table 15 shows that for homogeneous 

6 It may be tempting to provide an intuitive interpretation of this conclusion by trying to relate this lack of change of £ with 
the term ] ^_ de ^' r ^ v t ) of Equation (1) that bounds the total shock propagated between successive time steps. This explanation 
may not be fully correct for at least two reasons. Firstly, the above conclusion does not hold for homogeneous networks 
which also use Equation (1) for insolvency cascading. Secondly, in our simulations we have found this term to limit insolvency 
cascading only in a minority of cases. 
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in-arborescence networks the vulnerability index under coordinated shocks is at least as much as the vulner- 
ability index under idiosyncratic shocks 84.62% of the time, and for (0.1, 0.95)-heterogeneous ER networks 
of average degree 6 the vulnerability index under coordinated shocks is at least as much as the vulnerability 
index under idiosyncratic shocks 98.99% of the time 7 . Thus, we conclude: 

(F) correlated shocking mechanisms are more appropriate to measure the worst- case 
stability compared to idiosyncratic shocking mechanisms. 

4.5 Phase transition properties of stability 

Phase transitions are quite common when one studies various topological properties of graphs [43]. The 
vulnerability index £ shows several sharp phase transitions 8 . Two such interesting cases are discussed below. 

4.5.1 Dense homogeneous networks 

Based on the behavior of £ with respect to $ — 7, we observe that, for smaller value of k and for denser ER 
and SF networks under coordinated or idiosyncratic shocks, there is often a sharp decrease of stability when 
7 was decreased beyond a particular threshold value. For example, with $ = 0.5 and K. = 0.1, both 100 
node dense (average degree 6) SF and ER homogeneous networks exhibited more than ninefold increase in £ 
around 7 = 0.15 and 7 = 0.1, respectively 9 . 

To investigate the extent of such a sharp decrease around a threshold value of 7 in the range [0.05, 0.2] , 
we computed, for each of the five homogeneous network types under coordinated shocks and for each values 
of the parameters |V|, $, E /i and k, , the ratio 

max {£) — min {£) 

. , _ , 0.05<7<0.2 1 ' 0.05<7<0.2 1 ' 

A (n, $, E /i, «) = — m ^ m 

max (£j — mm (£j 

entire range of 7 entire range of 7 

that provides the maximum percentage of the total change of the vulnerability index that occurred within 
this narrow range of 7. The values of A (n, <£>, E /i, k) for the dense ER and SF homogeneous networks under 
coordinated shocks are shown in supplemental Table 16 for $ = 0.5,0.8 and k = 0.1,0.2,0.3,0.4,0.5 (the 
behaviour of £ is similar for other intermediate values of $). If the growth of £ with respect to 7 was uniform 
(linear) or near uniform over the entire range of 7, A would be approximately A = q'^^q ^ = 0.375; thus, 
any value of A significantly higher than A indicates a sharp transition within the above-mentioned range of 
values of 7. As Table 16 shows, a significant majority of the entries for $ < 0.8 and k < 0.2 are 2 A or more. 

4.5.2 Homogeneous in-arborescence networks 

Homogeneous in-arborescence networks under coordinated shocks (and to a lesser extent under idiosyncratic 
shocks) exhibited a sharp increase in stability as the ratio E /i of the total external asset to the total interbank 
exposure the system is increased beyond a particular threshold provided the equity to asset ratio 7 was 
approximately 50% of the shock parameter $. For example 10 , for a 50-node homogeneous in-arborescence 
network under coordinated shock, £ exhibited a sharp decrease from 0.76 to 0.18 for E /i G [0.75, l] , tC = 0.1, 
$ = 0.5 and 7 = 0.25 = */2. 



7 Visual illustrations of Conclusion (F) are shown in supplemental Fig. S6- Sll. 
8 For visual illustratation, see supplemental Fig. SI- S4. 
9 See supplemental Fig. SI for a visual illustration. 
10 Sec supplemental Fig. S4 for a visual illustration. 



12 



To investigate the extent of such a sharp decrease of £ around a threshold value of E /i in the range [0.5, l] 
with 7 « */2, we computed, for each type of shocking mechanism, and for each values of the parameters n, 
$, 7 « */2, and /t of the homogeneous in-arborescence network, the ratio 

max {£) — min {£) 

. , _ , 0.5<E//<1 1 ' 0.5<E//<1 1 J 

A(n,$,7,/s) = - 



max {£} — min {£} 

entire range of E entire range of E 



that provides the maximum percentage of the total change of the vulnerability index that occurred within 
this range of E /i. If the growth of £ with respect to E /i was uniform (linear) or near uniform over the entire 
range of E /i 1 A would be approximately S — 3 \ Z° 5 2 5 ~ 0.16; thus, any value of A significantly higher than 
S indicates a sharp transition within the above-mentioned range of E /i. As Table 17 shows, when 7 = */2 
a significant majority of the entries are coordinated shocks and many entries under idiosyncratic shocks are 
at least 2 S. 

An informal intuition behind the type of sharp decrease of £ is as follows. Suppose that at time i = 0a 
node u is shocked and becomes insolvent. For an in-arborescence network, the amount of shock transmitted 

from u to each of its in- neighbors is a « ($ — 7) x [ 1 H ■ — ) . Suppose that an in- neighbor v of 

\ n in-degree(u) / 

u has in-degree(w) = 0. Then, c v (0) = 7 x (E/n) and 

(E \ E 

1 + — -1 TT -7X - 
n m-degree(u) / n 

Assuming 7 w */2, we have 

/EE 
i/«7X 1 + 



n in-degree(u) n 

Thus, assuming in-degree(u) > 1, beyond the threshold value of E = (l — - n degree(tt) ) x n -< v wou ld be 

strictly negative, the node v will not become insolvent at any time t > 1, and the shock will not propagate 
any further through v. If u has many such in-neighbors of zero in-degrees, then the stability of the network 
will improve sharply around such a threshold value of E. If the in-degrees of the in-neighbors of u are not 
zero but bounded by a small constant, the same effect will be produced around a different threshold value 
oiE. 



5 Concluding remarks 

In this paper, we have initiated a methodology for systematic investigation of the global stabilities of financial 
networks that arise in OTC derivatives market and elsewhere. Our results can be viweed as a much needed 
beginning of a systematic investigation of these issues, with future research works concentrating on further 
improving the network model, the stability measure and parameter choices. 
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P = 0.6 


/3 = 0.95 P = 0.6 


P = 0.95 


P = 0.6 


/3 = 0.95 (9 = 0.6 


coordinated shock 


66.91% 60.22% 


99.26% 


98.91% 


98.46% 98% 


98.22% 


91.68% 


99.13% 97.4% 


idiosyncratic shock 


92.75% 81.79% 


97.76% 


96.81% 


98.16% 97.61% 


98.86% 


94.84% 


98.83% 97.22% 



Table 1: Comparison of stabilities of (a, /3)-heterogeneous networks with their homogeneous counter-parts 
over all parameter ranges. The numbers are the percentages of data points for which & a ,i3) -heterogeneous was 

at least ^homogeneous- 



average values of | max .25<E/r<3.5 {£} - min . 25 <E/i<3.5 {£} | 



under coordinated shock 


under idiosyncratic shock 


(0.1, 0.95)-heterogeneous in-arborescence 


0.017 


0.045 


(0.2, 0.6)-heterogeneous in-arborescence 


0.007 


0.017 


(0.1, 0.95)-heterogeneous ER, average degree 3 


0.066 


0.073 


(0.2, 0.6)-heterogeneous ER, average degree 3 


0.040 


0.041 


(0.1, 0.95)-heterogeneous ER, average degree 6 


0.111 


0.116 


(0.2, 0.6) -heterogeneous ER, average degree 6 


0.084 


0.078 


(0.1, 0.95)-heterogeneous SF, average degree 3 


0.119 


0.094 


(0.2, 0.6)-heterogeneous SF, average degree 3 


0.034 


0.032 


(0.1, 0.95)-heterogeneous SF, average degree 6 


0.200 


0.179 


(0.2, 0.6)-heterogeneous SF, average degree 6 


0.054 


0.054 



Table 2: Absolute values of the largest change of the vulnerability index £ in the range 0.25 < E /i < 3.5. 
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t = ; V; nS oivont(0) = ; continue =TRUE 

for every node v, let in-degree(«, 0) = in-degree(u) 

let V' be the set of nodes selected for initial insolvency 

for every node u £ V do 

c u (0) = c u — $e u 
endfor 



while ( ( continue =TRUE) A ( insolvent (t) + V)) do 
for every node u£F \ insolvent (t) do 
c u (t + 1) = c u (t) 



for every node w £ V \ V| n 



t (t) do 



if 



(cv(t) < 0) A ((«.«) G F )) then 

{ | c (*) | , b„ } 



in-degree(t>, i) 



decrease c u (£ + 1) by — 

endif 
endfor 

Vjnsolvcnt(^ ~t~ 1) — insolvent (^) 



( if c u (i) < then 
< add the node u to Insolvent (i + 1) 
y endif 
endfor 

if ( insolvent (t + 1) = V inso lvont (t ) ) then 

continue = FALSE 
endif 



rooti 




Figure 2: An in-arborescence graph. 



increase t by 1 
for every node u £V \ ^insolvent (t) do 

in-degree(it, t) = ((«,«) £F)A(d G V\ insolvent (t))}| 

endfor 



endwhile 



Figure 1: Pseudo-code of insolvency propagation method. 
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(0,1,0.95) ER 


(0.2,0.6) ER 


(0,1,0.95) SF 


(0.2,0.6) SF 


average degree 6 


average degree 6 


average degree 6 


average degree 6 


versus 


versus 


versus 


versus 


average degree 3 


average degree 3 


average degree 3 


average degree 3 


coordinated idiosyncratic 


coordinated idiosyncratic 


coordinated idiosyncratic 


coordinated idiosyncratic 


shock shock 


shock shock 


shock shock 


shock shock 


89.3% 82.39% 


68.12% 61.46% 


85.51% 73.81% 


69.29% 73.07% 



(A) 



(0,1,0.95) SF average degree 3 and average degree 6 
versus 

(0.1,0.95)-heterogeneous in-arborescence (SF average degree 1) 



(0.2,0.6) SF average degree 3 and average degree 6 
versus 

(0.2,0.6)-heterogeneous in-arborescence (SF average degree 1) 



coordinated shock 
85.7% 



idiosyncratic shock 
81.86% 



coordinated shock 
56.21% 



idiosyncratic shock 
51.07% 



(B) 



Table 3: Effect of connectivity on the global stability under coordinated and idiosyncratic shocks for 
(A) (a, /3)-heterogeneous ER and SF networks and (B) (a, /3)-heterogeneous in-arborescence versus (a, /3)- 
heterogeneous SF networks. The number (percentage) shown for a comparison of the type "network A versus 
network B" indicates the percentage of data points for which the stability of networks of type A was at least 
as much as that of networks of type B. 
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Table 4: Residual instabilities of homogeneous versus heterogeneous networks under coordinated shocks. 
The percentages shown are the percentages of networks for which £ < 0.05 or £ < 0.1 or £ < 0.2. 
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ER, 


average degree 3 


66% 


100% 


100% 


26% 


74% 


100% 






homogeneous 


ER, 


average degree 6 


100% 


100% 


100% 


100% 


100% 


100% 








SF, 


average degree 3 


29% 


72% 


100% 


19% 


53% 


88% 








SF, 


average degree 6 


100% 


100% 


100% 


88% 


100% 


100% 










in-arborescence 


0% 


1% 


12% 


0% 


1% 


11% 








ER, 


average degree 3 


0% 


0% 


15% 


0% 


1% 


10% 


\v\ 


= 100 


(0.1, 0.95)-heterogeneous 


ER, 


average degree 6 


6% 


7% 


16% 


6% 


6% 


10% 








SF, 


average degree 3 


0% 


6% 


23% 


0% 


0% 


16% 








SF, 


average degree 6 


8% 


16% 


30% 


6% 


12% 


19% 










in-arborescence 


0% 


1% 


11% 


0% 


1% 


11% 








ER, 


average degree 3 


6% 


10% 


18% 


0% 


7% 


16% 






(0.2, 0.6)-heterogeneous 


ER, 


average degree 6 


7% 


12% 


22% 


6% 


8% 


17% 








SF, 


average degree 3 


0% 


2% 


18% 


0% 


0% 


16% 








SF, 


average degree 6 


5% 


9% 


18% 


2% 


6% 


16% 










in-arborescence 


76% 


78% 


81% 


36% 


67% 


81% 








ER, 


average degree 3 


76% 


100% 


100% 


22% 


73% 


93% 






homogeneous 


ER, 


average degree 6 


100% 


100% 


100% 


100% 


100% 


100% 








SF, 


average degree 3 


26% 


52% 


100% 


20% 


42% 


88% 








SF, 


average degree 6 


100% 


100% 


100% 


87% 


100% 


100% 










in-arborescence 


0% 


1% 


12% 


0% 


1% 


12% 








ER, 


average degree 3 


0% 


2% 


16% 


0% 


0% 


5% 


\v\ 


= 300 


(0.1, 0.95)-heterogeneous 


ER, 


average degree 6 


0% 


6% 


12% 


0% 


2% 


6% 








SF, 


average degree 3 


0% 


9% 


24% 


0% 


1% 


16% 








SF, 


average degree 6 


7% 


16% 


28% 


6% 


9% 


19% 










in-arborescence 


0% 


1% 


11% 


0% 


1% 


11% 








ER, 


average degree 3 


6% 


8% 


16% 


0% 


2% 


16% 






(0.2, 0.6)-heterogeneous 


ER, 


average degree 6 


6% 


8% 


17% 


6% 


7% 


16% 








SF, 


average degree 3 


0% 


1% 


16% 


0% 


0% 


16% 








SF, 


average degree 6 


2% 


6% 


16% 


2% 


5% 


16% 



Table 5: Residual instabilities of homogeneous versus heterogeneous networks under idiosyncratic shocks. 
The percentages shown are the percentages of networks for which £ < 0.05 or £ < 0.1 or £ < 0.2. 
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coordinated shock 















0.6,7 = 


55 


$ = 


0.6,7 = 


150 












i < 0.05 


£ < 0.1 


£ < 0.2 


£ < 0.05 


£ < 0.1 


£ < 0.2 










in-arborescence 


93% 


93% 


93% 


2% 


69% 


93% 








ER, 


average degree 3 


97% 


100% 


100% 


64% 


95% 


100% 






homogeneous 


ER, 


average degree 6 


100% 


100% 


100% 


100% 


100% 


100% 








SF, 


average degree 3 


56% 


96% 


100% 


44% 


87% 


100% 








SF, 


average degree 6 


100% 


100% 


100% 


100% 


100% 


100% 










in-arborescence 


0% 


0% 


1% 


0% 


0% 


0% 








ER, 


average degree 3 


0% 


0% 


5% 


0% 


0% 


1% 


\v\ = 


50 


(0.1, 0.95)-heterogeneous 


ER, 


average degree 6 


9% 


11% 


14% 


5% 


7% 


7% 








SF, 


average degree 3 


5% 


9% 


19% 


0% 


3% 


8% 








SF, 


average degree 6 


22% 


28% 


41% 


10% 


15% 


17% 










in-arborescence 


0% 


0% 


9% 


0% 


0% 


9% 








ER, 


average degree 3 


5% 


7% 


20% 


3% 


7% 


16% 






(0.2, 0.6)-hcterogeneous 


ER, 


average degree 6 

o o 


11% 


14% 


26% 


7% 


7% 


19% 








SF, 


average degree 3 


2% 


8% 


23% 


1% 


3% 


18% 








SF, 


average degree 6 


9% 


15% 


26% 


7% 


8% 


18% 










in-arborescence 


93% 


93% 


93% 


23% 


84% 


93% 








ER, 


average degree 3 


85% 


100% 


100% 


47% 


86% 


100% 






homogeneous 


ER, 


average degree 6 


100% 


100% 


100% 


100% 


100% 


100% 








SF, 


average degree 3 


37% 


75% 


100% 


35% 


71% 


100% 








SF, 


average degree 6 


100% 


100% 


100% 


96% 


100% 


100% 










in-arborescence 


0% 


0% 


0% 


0% 


0% 


0% 








ER, 


average degree 3 


0% 


1% 


7% 


0% 


0% 


3% 


\v\ = 


100 


(0.1, 0.95)-heterogeneous 


ER, 


average degree 6 


7% 


7% 


9% 


7% 


7% 


7% 








SF, 


average degree 3 


0% 


1% 


9% 


0% 


0% 


5% 








SF, 


average degree 6 


10% 


15% 


21% 


7% 


7% 


11% 










in-arborescence 


0% 


0% 


9% 


0% 


0% 


9% 








ER, 


average degree 3 


0% 


7% 


17% 


0% 


5% 


16% 






(0.2, 0.6)-heterogcncous 


ER, 


average degree 6 


7% 


7% 


19% 


7% 


7% 


16% 








SF, 


average degree 3 


1% 


3% 


15% 


0% 


1% 


13% 








SF, 


average degree 6 


7% 


9% 


18% 


7% 


7% 


17% 










in-arborescence 


93% 


93% 


93% 


79% 


93% 


93% 








ER, 


average degree 3 


93% 


100% 


100% 


48% 


85% 


100% 






homogeneous 


ER, 


average degree 6 


100% 


100% 


100% 


100% 


100% 


100% 








SF, 


average degree 3 

o o 


28% 


56% 


97% 


28% 


56% 


97% 








SF, 


average degree 6 


100% 


100% 


100% 


93% 


100% 


100% 










in-arborescence 


0% 


0% 


0% 


0% 


0% 


0% 








ER, 


average degree 3 


0% 


1% 


7% 


0% 


0% 


3% 


\v\ = 


300 


(0.1, 0.95)-heterogeneous 


ER, 


average degree 6 


7% 


7% 


7% 


3% 


7% 


7% 








SF, 


average degree 3 


0% 


0% 


11% 


0% 


0% 


3% 








SF, 


average degree 6 


7% 


7% 


18% 


7% 


7% 


7% 










in-arborescence 


0% 


0% 


9% 


0% 


0% 


9% 








ER, 


average degree 3 


1% 


7% 


17% 


1% 


5% 


16% 






(0.2, 0.6)-heterogeneous 


ER, 


average degree 6 


7% 


7% 


16% 


7% 


7% 


16% 








SF, 


average degree 3 


0% 


1% 


13% 


0% 


0% 


13% 








SF, 


average degree 6 


7% 


7% 


16% 


7% 


7% 


16% 



Table 6: Residual instabilities of homogeneous versus heterogeneous networks under coordinated shocks. 
The percentages shown are the percentages of networks for which £ < 0.05 or £ < 0.1 or £ < 0.2. 
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coordinated shock 















n 7 ~ — 
u. 1 , 7 — 


J.UO 


* — 


u. 1, 7 — 


J.UU 












f <r fl D5 




<^ u.^ 
















in-arborescence 


if y* /C 




QS% 


00 /(J 


88% 


QS% 








ER, 


average degree 3 


97% 


100% 


100% 


81% 


100% 


100% 






I1U II lUgCIlCU Lib 


ER, average degree 6 


1 00% 

1UU /o 


1 00% 


1 00% 

1UU /o 


1 00% 

J- LFLF /0 


1 00% 


1 00% 








SF, 


average degree 3 


5fi% 


Qfi% 


1 00% 

luu /(; 


^2% 


QR% 


1 00% 

1 V M J /(j 








SF, 


average degree 6 


J. WW /0 


1 00% 

-L WW /0 


1 00% 

J- VJ VJ /o 


1 00% 


1 00% 


1 00% 

J- WW /o 










in-arborescence 


u /o 


1 /O 


7% 
f /o 


u /o 


u/o 


u /o 








ER, average degree 3 


0% 


0% 


5% 


0% 


0% 


1% 
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\ v \ - 
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- o\j 
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I U.l, U.i?U LLT UgLllLW Lib 
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±u /o 


13% 

-LO /0 


1fi% 

-i-U /u 


7% 


7% 


Q% 








SF, 


average degree 3 


7% 
i /o 


19% 


25% 


1 % 


4% 


9% 








SF, 


average degree 6 




39% 
/o 


4Q% 


14% 


1 5% 

J-U /o 


21% 

^ -L /0 










in-arborescence 


0% 

u /o 


0% 

u /o 


10% 


0% 

U/O 


0% 

u /o 


9% 

/0 








ER, 


average degree 3 


5% 


8% 


21% 


3% 


7% 


17% 






I l-Viptprncypripniic; 

l U . Z; , U . U f LV^l U^LllLU Llo 


ER, average degree 6 


1 3% 


18% 

-LO /0 


30% 


7% 


7% 

1 /u 


20% 

^U /u 
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jivpt'Jio'p Hpctpp X 


3% 


9% 


25% 


1% 

J- /(J 


3% 

O /O 


1Q% 

-L /0 








SF 
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U /{J 


16% 

J-U /o 


28% 


8% 

O /(J 


9% 

/0 


19% 

J- cf /O 
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111 dl UU1LSLL11LL 


yo /o 


yo /o 


yo /o 


Ol/O 


yo /o 


yo /o 








ER, 


average degree 3 


85% 
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100% 


62% 


97% 


100% 






HWlIlUJi^llLW Lib 


FR 

1 1 \ . 


cxVL-Ielf^L Ll^^l U 


1 00% 

J- LFLF /0 


1 00% 

J- LFLF /0 


1 00% 


1 00% 

J- LFLF /0 


1 00% 


1 00% 








SF 


tii^pvft or 1 Hocfror 1 
dvl^ldtil^ til o 


S7% 


75% 

1 *J /(J 


1 00% 

_L VJ VJ /(J 


S7% 

O J /(J 


75% 


1 00% 

_L VJ VJ /O 








SF 


dVL-IelJ^L UC^ICC U 


1 00% 

-L WW /0 


1 00% 


1 00% 


1 00% 

-L LFLF /0 


1 nn% 

J- vvvv /o 


1 nn% 

J- WW /o 










1 y | ■! 1' 1 1/ lf( i^'/'iH 1 / ■ ( . 

111-ai uoiescenct 


n% 

u /o 


0% 
u /o 


u /o 


n% 

u /o 


u /o 


u /o 








ER, 


average degree 3 


0% 


1% 


7% 


0% 


0% 


5% 


\v\ - 

| V | — 


1 on 


1(1 1 (1 Qm n L pr pi ~n pAi i q 
I U.l, \J.jO J LLT UgL-llLW Lib 


ER, 


average degree 6 


7% 
I /o 


7% 
I /o 


Q% 


7% 


7% 


7% 








SF, 


average degree 3 


0% 

u /o 


2% 


1 5% 


0% 

U/O 


0% 
u /o 


5% 
u /o 








SF, 


average degree 6 


11% 

j. j. /u 


15% 

-L'-J /O 


27% 


7% 

1 /O 


7% 

1 /U 


11% 

-L -L /U 










in-arborescence 


0% 

u /o 


0% 

u /o 


q% 

C7 /(J 


0% 

u/o 


0% 
u /o 


9% 

J /0 








ER, 


average degree 3 


0% 


7% 


18% 


0% 


7% 


16% 






(0 9 i- np1"PT - no'PTi pmiQ 

I U . Zi, U . U 1 11L- LLI UgLllLU LIS 


ER, 


average degree 6 


7% 


7% 


19% 


7% 

1 /o 


7% 

( /O 


16% 

J- U /O 








SF, 


average degree 3 


1 /o 


4% 


17% 

J- 1 /V 


0% 

u/o 


2% 
^ /o 


1 3% 

1O/0 








SF, 


average degree 6 


/o 


10% 

-LU /0 


19% 


7% 

1 /o 


8% 
/o 


17% 

1 ( /o 










in-arborescence 






t/O /U 


OLF /O 


yo /o 


yo /o 








ER, 


average degree 3 


93% 


100% 


100% 


65% 


97% 


100% 






HWlIlUgLULW Lib 


ER, 


average degree 6 


1 00% 


1 00% 


1 00% 


1 00% 


1 00% 


1 00% 








SF, 


average degree 3 


28% 


56% 

OLF /(J 


Q7% 


2S% 


56% 


97% 

<J 1 /(J 








SF, 


average degree 6 


1 00% 


1 00% 


1 00% 
J- vvvv /o 


1 00% 

-L LFLF /O 


1 00% 

J- WW /o 


1 00% 

J- WW /o 










in-arborescence 


u/o 


u /o 


u /o 


u /o 


u/o 


u /o 








ER, 


average degree 3 


0% 


1% 


7% 


0% 


0% 


4% 


\v\ = 


300 


(0.1, 0.95)-heterogeneous 


ER, 


average degree 6 


7% 


7% 


7% 


6% 


7% 


7% 








SF, 


average degree 3 


0% 


0% 


15% 


0% 


0% 


3% 








SF, 


average degree 6 


7% 


7% 


21% 


7% 


7% 


7% 










in-arborescence 


0% 


0% 


9% 


0% 


0% 


9% 








ER, 


average degree 3 


1% 


7% 


17% 


0% 


7% 


16% 






(0.2, 0.6)-heterogeneous 


ER, 


average degree 6 


7% 


7% 


16% 


7% 


7% 


16% 








SF, 


average degree 3 


0% 


1% 


14% 


0% 


0% 


13% 








SF, 


average degree 6 


7% 


8% 


16% 


7% 


7% 


16% 



Table 7: Residual instabilities of homogeneous versus heterogeneous networks under coordinated shocks. 
The percentages shown are the percentages of networks for which £ < 0.05 or £ < 0.1 or £ < 0.2. 
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1% 


6% 
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19% 
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10% 

XU /0 


0% 

U/O 
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5% 
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7% 
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Table 8: Residual instabilities of homogeneous versus heterogeneous networks under coordinated shocks. 
The percentages shown are the percentages of networks for which £ < 0.05 or £ < 0.1 or £ < 0.2. 
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Table 9: Residual instabilities of homogeneous versus heterogeneous networks under coordinated shocks. 
The percentages shown are the percentages of networks for which £ < 0.05 or £ < 0.1 or £ < 0.2. 
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Table 10: Residual instabilities of homogeneous versus heterogeneous networks under coordinated shocks. 
The percentages shown are the percentages of networks for which £ < 0.05 or £ < 0.1 or £ < 0.2. 
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Table 11: Residual instabilities of homogeneous versus heterogeneous networks under coordinated shocks. 
The percentages shown are the percentages of networks for which £ < 0.05 or £ < 0.1 or £ < 0.2. 
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Table 12: Residual instabilities of homogeneous versus heterogeneous networks under coordinated shocks. 
The percentages shown are the percentages of networks for which £ < 0.05 or £ < 0.1 or £ < 0.2. 
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Table 13: Residual instabilities of homogeneous versus heterogeneous networks under coordinated shocks. 
The percentages shown are the percentages of networks for which £ < 0.05 or £ < 0.1 or £ < 0.2. 
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ER average degree 3 


versus average degree 6 


SF average degree 3 versus average degree 6 


coordinated shock 


idiosyncratic shock 


coordinated shock idiosyncratic shock 


97.43% 


97.05% 


98.89% 98.29% 



Table 14: Effect of connectivity on the global stability for homogeneous networks under coordinated and 
idiosyncratic shocks. The number (percentage) shown for a comparison of the type "network A versus 
network B" indicates the percentage of data points for which the stability of networks of type A was at least 
as much as that of networks of type B. 



In-arborescence 
a = 0.1 a = 0.2 
P = 0.95 = 0.6 

56.64% 57.27% 



(a, /3)-heterogeneous networks 



ER average degree 3 
a = 0.1 a = 0.2 
P = 0.95 P = 0.6 

89.66% 90.97% 



ER average degree 6 
a = 0.1 a = 0.2 
P = 0.95 P = 0.6 

98.99% 98.04% 



SF average degree 3 
a = 0.1 a = 0.2 
P = 0.95 P = 0.6 
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SF average degree 6 
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homogeneous networks 
ER average degree 6 
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SF average degree 3 
81.15% 



SF average degree 6 
54.80% 



Table 15: Comparisons of strengths of coordinated versus idiosyncratic shocks. The numbers (percentages) 
indicate the percentage of total number of data points (combinations of parameters 7, E and k) for 
that network type that resulted in £ c being at least £ r , where £ c and £ r are the vulnerability indices under 
coordinated and idiosyncratic shocks, respectively. 
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0.66 
0.63 
0.6 



0.59 
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0.56 



0.78 
0.79 
0.87 



0.79 
0.81 
0.8 



0.79 

0.8 

0.78 



0.76 
0.77 
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0.72 
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0.71 



0.98 
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0.94 
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0.84 
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0.72 



0.77 

0.68 
0.61 
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0.63 
0.58 
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0.78 
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0.79 



0.76 

0.74 
0.75 



0.69 
0.66 
0.68 



0.59 
0.57 
0.6 



0.94 
0.92 
0.77 



0.87 
0.77 

0.62 



0.76 

0.69 
0.61 



0.66 
0.63 
0.6 



0.58 
0.56 
0.55 



0.78 
0.79 
0.87 
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0.81 
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0.77 
0.75 



0.72 
0.73 
0.71 



0.98 
0.98 
0.94 



0.94 
0.87 
0.84 



0.86 

0.73 
0.72 



0.77 

0.68 
0.61 



0.7 
0.63 
0.58 



0.8 

0.81 

0.81 



0.79 
0.79 
0.79 



0.76 

0.74 
0.74 



0.68 
0.65 
0.67 



0.57 
0.55 
0.58 



0.94 
0.92 
0.77 



0.88 
0.77 

0.61 



0.75 

0.69 
0.61 



0.65 
0.62 
0.59 



0.55 
0.53 
0.52 



0.78 
0.79 
0.87 



0.79 
0.81 
0.8 



0.79 

0.8 

0.78 



0.76 
0.77 
0.75 



0.72 
0.73 
0.71 



0.98 
0.98 
0.94 



0.94 
0.87 
0.84 



0.86 

0.73 
0.72 



0.77 

0.68 
0.61 



0.7 
0.63 
0.58 



0.81 
0.81 
0.81 



0.78 
0.79 
0.78 



0.75 

0.73 
0.73 



0.67 
0.64 
0.65 



0.55 
0.54 
0.56 



0.94 
0.92 
0.76 



0.87 
0.76 

0.61 



0.75 

0.69 
0.61 



0.64 
0.62 
0.59 



0.54 
0.52 
0.52 



max {£} — min {£} 

Table 16: Values of A (n, $, E/I, k) — 05 <7<o.2 -, , °'° 5 ~ 7 ~°' 2 for homogeneous dense ER and 

max |£j — min |£j 

entire range of 7 entire range of 7 

SF networks under coordinated shocks. Entries that are at least 2 x 0.375 arc shown in boldface. 
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max {£} — min {£} 

Table 17: Values of A(n, $, 7, k) = ^ = — ^ — ^ j^j- for homogeneous in-arborescence 

max u] - min |£j 

entire range of E entire range of E 

networks. Entries that are at least 2 x 0.16 are shown in boldface black. Entries that are at least | x 0.16 
are shown in boldface gray. 
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APPENDIX 




A Illustration of calculations of parameters for a sample financial 
network shown in Fig. 3 

n = number of nodes = 5 
m = number of edges = 7 
7 = 0.1 

I = total inter-bank exposure = m = 7 
E = total external asset = 21 = 14 

Figure 3: A banking network. 

(a) Homogeneous version of the network 

• w = weight of every edge = I /m = 1. 

• E — E — E — E — E — 14 A 

• ^V2 — — 1; ^V3 ^. 

• b Vl = in-degree(ui) = 2,b V2 = in-degree(w 2 ) = l,b V3 = in-degree(u 3 ) = l,b V4 = in-degree(u 4 ) = 
3,b V5 = in-degree^) = 0. 



-■V 5 — \"V 5 ^V 5 J ~ J-'Vz 



• e Vl = 




(K- 


i Vi ) + E V4 


• a Vl = 


b Vl + E Vl 



*V 3 — "t> 3 ' J-/V3 



l V 5 



b V5 + E V5 — 2.8. 

• c Vl =^a Vl = 0.48, c V2 ='ja V2 = 0.38, c V3 = ja V3 = 0.38, c„ 4 = -ya V4 = 0.58, c„ 5 = ja V5 = 0.28. 
(b) Heterogeneous version of the network 

Suppose that 95% of E is distributed (equally) on the two banks vi and w 2 , and the rest 5% of E is distributed 
(equally) on the remaining three banks. Thus: 



3 = 1 



= 6.65, E V2 = a V2 \ E - ipvj ~ Lv i) 

i=i 



6.65, E V3 — a V3 \ E - - S 

3 = 1 



E V4 = a V4 y E ~JL - S) J ~ °- 233 ' E v* = <*V5 \^~J2 ( bv * ~ ^ ) J ~ °- 233 

Suppose that 95% of / is distributed (equally) on the three edges f\ = («2, «i), fi = W4), f% = («4, W2), and 
the remaining 5% of / is distributed (equally) on the remaining four edges f^ — («3,vi),/5 = («3, V4), ,/6 = 
{v 5 ,v 4 ),f7 = (v 5 ,v 3 ). Then, 

w(A) = w(f 2 ) = w(h) = » 2.216, w(U) = w(f 5 ) = w(f 6 ) = w(f 7 ) = = 0.08725 



32 



, , . b Vl = wifi) + w(f 4 ) m 2.30325, t Vl = w{f 2 ) =2.216, e Vl = {b Vl - l Vi ) + E Vl ps 6.7365 
for bank V\: 

dvi = b vi + E V1 = 8.9525, c vi = ja vi = 0.8925 

b V2 = w(f 3 ) w 2.216, l V2 = w(fi) « 2.216, e„ 2 = (6„ 2 - + E V2 = 6.65 
tor bank v 2 : ' 
a>v 2 = b V2 + ~ 8.866, c„ 2 = 7«„ 2 w 0.8666 

6 TO = H/t) = 0.08725, V 3 = w{U) + w(f 5 ) = 0.1745, e V3 = (b V3 - i V3 ) + E V3 w 0.14575 

tor bank w 3 : 

a„ 3 = b V3 + E V3 « 0.32025, c„ 3 = 7a„ 3 w 0.032035 
6„ 4 = w(/ 2 ) + w(/ 5 ) + w(f 6 ) w 2.39050, ^ = w(/ 3 ) » 2.216, e„ 4 = (6 W - v) + ^ w » 0.4075 

for bank V4 : 

a Vi = b Vi + E V4 » 2.6235, c„ 4 = 7 a„ 4 w 0.26235 

, , . &„ 5 =0, i V5 =w{f 6 )+w{f 7 ) = 0.1745, e„ 5 = (&„„ -^J+E^ ss 0.0585 
for bank v 5 : 

a V5 = b Vb + E Vb ps 0.233, c V5 = 7 a„ 5 ps 0.0233 
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Supplementary color figures SI — Sll 
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ER model (average degree 6), coordinated shock 



ER model (average degree 3), coordinated shock 



SF model (average degree fi), coordinated shock SF model (average degree 3), coordinated shock ^^^m 

in-arborescence (average degree 1), coordinated shock ■■■ 



E = \V\ = 50 K = 0.1 * = 0.5 E = \V\ = 50 JC = 0.5 * = 0.5 E = |V"| = 50 K = 0.9 * = 0.5 




-y 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 * - -y 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 - -y 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 

7 0.45 0.40 0.35 0.30 0.25 0.20 0.15 0.10 0.05 7 0.45 0.40 0.35 0.30 0.25 0.20 0.15 0.10 0.05 7 0.45 0.40 0.35 0.30 0.25 0.20 0.15 0.10 0.05 



(A) (B) (C) 



E=\V\ = 100 fC = 0.1 * = 0.5 E = \V\ = 100 K = 0.5 * = 0.5 E = \V\ = 100 K = 0.9 * = 0.5 




$ — 7 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 * - -7 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 * — 7 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 

7 0.45 0.40 0.35 0.30 0.25 0.20 0.15 0.10 0.05 7 0.45 0.40 0.35 0.30 0.25 0.20 0.15 0.10 0.05 7 0.45 0.40 0.35 0.30 0.25 0.20 0.15 0.10 0.05 
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Figure SI: (drawn in color) Effect of variations of equity to asset ratio (with respect to shock) on the 
vulnerability index £ for homogeneous networks. Lower values of £ imply higher global stability of a network. 
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Figure S2: (drawn in color) Effect of variations of equity to asset ratio (with respect to shock) on the 
vulnerability index £ for (a, /3)-heterogeneous networks. Lower values of £ imply higher global stability of a 
network. 
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Figure S3: (drawn in color) Effect of variations of equity to asset ratio (with respect to shock) on the 
vulnerability index £ for (a, /3)-heterogeneous networks. Lower values of £ imply higher global stability of a 
network. 
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Figure S4: (drawn in color) Effect of variations of the total external to internal asset ratio E/I on the 
vulnerability index £ for homogeneous networks. Lower values of £ imply higher global stability of a network. 
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Figure S5: (drawn in color) Effect of variations of the total external to internal asset ratio E /i on the 
vulnerability index £ for (a, /3)-heterogeneous networks. Lower values of £ imply higher global stability of a 
network. 
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Figure S6: (drawn in color) Effect of variations of equity to asset ratio (with respect to shock) on the 
vulnerability index £ for homogeneous networks. Lower values of £ imply higher global stability of a network. 
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Figure S7: (drawn in color) Effect of variations of equity to asset ratio (with respect to shock) on the 
vulnerability index £ for (a, /3)-heterogeneous networks. Lower values of £ imply higher global stability of a 
network. 
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Figure S8: (drawn in color) Effect of variations of equity to asset ratio (with respect to shock) on the 
vulnerability index £ for (a, /3)-heterogeneous networks. Lower values of £ imply higher global stability of a 
network. 
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Figure S9: (drawn in color) Effect of variations of the total external to internal asset ratio E /i on the 
vulnerability index £ for homogeneous networks. Lower values of £ imply higher global stability of a network. 
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Figure S10: (drawn in color) Effect of variations of the total external to internal asset ratio E /i on the 
vulnerability index £ for (a, /3)-heterogeneous networks. Lower values of £ imply higher global stability of a 
network. 
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Figure Sll: (drawn in color) Effect of variations of the total external to internal asset ratio E /i on the 
vulnerability index £ for (a, /3)-heterogeneous networks. Lower values of £ imply higher global stability of a 
network. 
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